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ABSTRACT 

The  theory  of  the  sound  pressure  scattered  from  a  fluid-filled  spherical  shell  immersed 
in  a  second  fluid  is  developed.  An  arbitrary  combination  of  fluids  may  be  used  to  fill 
the  shell,  and  both  the  sound  source  and  the  receiver  may  be  in  the  near  field.  The 
theory  is  used  to  calculate  the  pressure  field  in  both  the  interior  and  exterior  fluids  and 
is  extended  to  determine  the  time-dependence  of  the  reflected  signal  and  the  target 
strength.  The  dependence  of  these  quantities  as  a  function  of  frequency  and  pulse 
type;  shell  material,  diameter  and  wall  thickness;  fluid  combination;  temperature;  and 
distance  of  the  receiver  from  the  sphere  is  examined,  and  is  shown  to  be  sensitive  to 
most  of  these  parameters.  The  performance  of  spheres  filled  with  a  Freon-113™  - 
ethanol  or  Fluorinert  FC-72™  -  n-hexane  mixture  is  considered. 
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Calculations  of  the  Sound  Scattering  from  Fluid- 
Filled  Spherical  Shell  Sonar  Targets 


Executive  Summary 

The  theory  of  the  sound  pressure  scattered  from  a  fluid-filled  spherical  shell 
immersed  in  a  second  fluid  is  developed.  The  prime  reason  for  developing  the  theory 
is  to  aid  in  the  design  of  fluid-filled  spheres  of  known  target  strength,  the  target 
strength  being  a  measure  of  the  reflected  intensity.  An  arbitrary  combination  of 
fluids  may  be  used  to  fill  the  shell,  and  both  the  sound  source  and  the  receiver  may 
be  placed  near  to  or  far  from  the  sphere.  The  theory  is  used  to  calculate  the  pressure 
field  in  both  the  interior  and  exterior  fluids  and  is  extended  to  determine  the  time- 
dependence  and  intensity  of  the  reflected  signal. 

Not  only  is  the  target  strength  found  to  be  sensitive  to  a  number  of  material 
parameters,  but  with  a  given  selection  of  materials,  it  is  also  sensitive  to  both  the 
form  of  the  incoming  sound  waves  and  the  environment.  The  target  strength  varies 
greatly  for  small  changes  in  the  frequency  if  continuous  waves  are  used. 
Fortunately,  these  variations  can  be  largely  eliminated  by  the  use  of  pulsed  signals, 
which  effectively  average  over  many  frequencies. 

A  shell  of  aluminium  or  pyrex  glass  optimises  the  target  strength,  but  as  these  are 
soft  or  brittle  materials,  stainless  steel  might  be  preferable  as  it  is  more  robust.  As 
might  be  expected,  increasing  the  diameter  of  the  sphere  will  give  a  larger  target 
strength.  A  thin  shell  wall  will  also  lead  to  a  large  target  strength,  but  may  be  too 
difficult  to  manufacture  and  will  be  easily  damaged.  A  thickness  of  0.8  mm  might  be 
classed  as  a  reasonable  compromise.  The  target  strength  is  quite  sensitive  to  the 
thickness,  so  variations  in  the  wall  thickness  need  to  be  avoided  in  the 
manufacturing. 

Choice  of  the  optimum  combination  of  fluids  with  which  to  fill  the  shell  depends 
on  the  temperature.  Variations  of  a  few  degrees  will  make  only  a  modest  change  to 
the  target  strength,  but  a  combination  chosen  for  warm  waters  will  behave  quite 
differently  in  cool  waters.  For  frequencies  below  200  kHz  a  mixture  of  Freon-113™ 
and  ethanol,  or  Fluorinert  FC-72™  and  n-hexane  will  work  quite  well  if  the  ratio  of 
liquids  is  chosen  to  give  a  refractive  index  between  1.8  and  1.9.  At  frequencies  above 
250  kHz  the  same  liquids  mixed  to  give  a  refractive  index  between  1.6  and  1.7  will 
produce  high  target  strengths. 
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1.  Introduction 

For  the  calibration  of  sonar  systems,  spherical  targets  are  advantageous  in  that  they  are 
aspect  independent.  That  is,  they  scatter  the  same  intensity  of  sound  into  a  particular 
direction,  regardless  of  the  target's  orientation.  However,  a  major  limitation  is  that  a 
solid  sphere  or  air-filled  spherical  shell  is  usually  a  weak  scatterer  in  comparison  to 
non-aspect  independent  reflectors  made  up  of  multiple  corner  reflectors.  An  attractive 
alternative,  providing  both  aspect  independence  and  strong  scattering,  is  a  fluid-filled 
spherical  shell  in  which  the  filling  fluid  is  chosen  so  as  to  cause  focusing  of  the 
reflected  sound  in  the  backscatter  direction.  This  direction  is  important  experimentally. 

This  study  considers  such  sonar  targets  and  the  effect  the  various  parameters  have  on 
their  performance.  The  aim  is  to  enable  selection  of  appropriate  materials  and  to  show 
how  the  performance  of  the  target  alters  with  small  changes  in  the  material  or 
environmental  parameters. 

The  scattering  of  sound  from  spherical  objects  has  been  the  subject  of  some 
investigation  in  the  past.  Rayleigh  (1945)  considered  the  limiting  case  where  the 
scatterer  was  small  compared  to  the  wavelength.  Morse  (1936)  derived  a  solution  for 
rigid  immovable  spheres,  which  were  not  necessarily  small  compared  to  the 
wavelength.  Faran  (1951)  allowed  penetration  of  sound  into  solid  spheres;  Junger 
(1952)  considered  the  scattering  by  thin  hollow  elastic  shells.  Goodman  and  Stern 
(1962)  also  considered  elastic  spherical  shells,  but  with  the  internal  and  external  fluids 
the  same.  Hickling  (1964)  analysed  evacuated  metal  spheres  in  water. 

Apart  from  restrictions  on  the  interior  fluid,  in  each  of  these  studies  it  was  assumed 
that  the  sphere  was  in  the  far  field  of  the  projector.  However,  the  theory  enabling  the 
scattered  sound  pressure  to  be  calculated  will  be  developed  in  the  following  section 
without  these  limitations,  by  solving  the  wave  equation  for  a  fluid-filled  spherical  shell 
immersed  in  another  fluid.  A  point  source  of  sound  waves  and  a  point  isotropic 
receiver  are  located  in  the  external  fluid.  The  former  may  be  arbitrarily  placed,  but  the 
receiver  must  be  no  further  away  from  the  centre  of  the  sphere  than  the  projector. 


2.  Theory 


2.1  Continuous  Waves 

Consider  a  fluid  medium  of  infinite  extent  whose  density  is  p,  and  which  supports 
longitudinal  sound  waves  of  speed  Cj.  In  the  medium  is  a  fluid-filled  spherical  shell, 
with  its  centre  located  at  the  origin  of  the  spherical  coordinate  system,  as  shown  in 
Fig.  1.  The  inner  and  outer  radii  of  the  shell  are  a  and  b,  respectively.  The  shell  has  a 
density  p2  and  supports  both  longitudinal  and  transverse  waves  of  velocity 
C2-  2,  respectively.  The  interior  fluid  has  a  density  Pj  and  a  longitudinal  sound  speed  of 
Cj.  Located  at  and  0  =  0  with  respect  to  the  origin  is  a  point  source  of 
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continuous  sinusoidal  sound  waves.  The  incident  pressure  a  distance  R  away  in  the 
external  fluid  is 


Pi=Po 


R 


(1) 


where  is  the  amplitude  of  the  pressure  wave  1  m  from  the  source, /is  the  frequency 
of  the  wave  and  j  is  the  wavenumber  in  the  fluid,  ,  is  related  to/  by 


(2) 


The  displacement  component  of  the  incident  pressure  can  be  written  as 


Pi 


=  Pn 


R 


(3) 


Receiver 


Figure  1:  Geometry  for  projector,  receiver  and  fluid-filled  sphere. 


When  the  sound  waves  strike  the  sphere,  some  of  the  sound  will  be  transmitted  into 
the  shell  and  the  filling  fluid,  and  some  will  be  scattered  back  into  the  exterior 
medium.  The  prime  focus  of  this  study  will  be  to  derive  the  scattered  pressure  at  the 
listening  point  of  r  =  and  0  =  0^.  This  can  be  found  by  satisfying  the  boundary 
conditions  at  the  surfaces  of  the  shell.  For  a  fluid  and  shell  behaving  in  a  linear 
fashion,  a  shell  of  zero  viscosity  and  with  no  cavitation  occurring,  these  are: 

i)  the  pressure  in  the  fluid  must  be  equal  to  the  normal  component  of  stress  in  the 
shell  at  the  interfaces. 
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ii)  the  normal  (radial)  component  of  displacement  of  the  fluid  must  be  equal  to  the 
normal  component  of  displacement  of  the  shell  at  the  interfaces,  and 

iii)  the  tangential  components  of  shearing  stress  must  vanish  at  the  surfaces  of  the 
shell. 

Symbolically, 


Pi  =-[rr] 

^r,l  “  ^r,2 

[r0]=[r(l)]  =  O 


P3  =-['•/■] 


al  r  =  b  3  =  2 


at  r  =  a 


[a-0]  =  [r(j)]  =  0 


where  the  /?,  and  are  the  pressures  in  the  external  and  internal  fluids,  with  the 
former  being  the  sum  of  the  incident  and  scattered  pressures 

Pi=Pi+Ps-  (5) 

The  are  the  radial  components  of  the  displacement  in  the  various  media  and  [r  r], 

[r  0]  and  [r  ({)]  are  the  radial  and  tangential  components  of  the  stress  in  the  shell. 

The  pressures  in  the  fluids  can  be  found  by  solving 

=  (6) 

in  which  i  is  1  or  3.  This  is  done  in  Appendix  A  using  the  spherical  coordinate  system, 
separating  out  the  time  and  position  variables  and  noting  that  P  has  no  ())  dependence, 
to  yield 

P  =  p{r,Q).T{t)  (7) 

where  in  the  external  and  internal  fluids  /?(r,0)  is 

Pi=XP/(cos0jA;h^(^2.,i';)  (8) 


P3=EPi(cos0jfiJ;(^^3rJ 


respectively,  and  j;  are  the  spherical  Hankel  and  Bessel  functions  of  order  /;  is 
the  Legendre  function  of  order  /. 
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It  will  be  more  convenient  to  use  different  summations  at  a  later  stage.  Using  a  series 
expansion  given  by  Morse  and  Feshbach  (1953,  p.  1574),  the  incident  pressure  in 
equation  (3)  can  be  written  as 

A  =  -in^uX(-l)'(2/+l)h;(fcL,i'p)j/(^L.A>/(cose  J  for  0<r^<r^  (10) 

(=0 

where  the  listening  point  is  at  r  =  and  0  =  0^  with  respect  to  the  origin,  and  R,  r^, 
and  0,  are  related  by  the  cosine  rule: 

=  r/  +  Vp  -  COS0,.  (H) 

Using  an  expression  of  Morse  and  Ingard  (1968,  p.  419)  it  can  be  shown  that  at  the 
listening  point  of  r  =  r,  and  9  =  0,,  the  scattered  pressure  will  be 

A  =  /^,it,,,S(-l)'(2/+l)C,h;(^,,rJh;(/:,,rJP,(cos0j  for  0  <  r,  <  (12) 

1=0 

where  the  C,  are  coefficients  to  be  determined.  Once  these  coefficients  have  been 
obtained,  the  scattered  pressure  will  be  known,  which  is  the  main  focus  of  this  study. 
The  pressure  in  the  external  fluid,  as  given  in  equation  (5),  is  found  by  summing 
equations  (10)  and  (12).  The  pressure  in  the  internal  fluid  is  of  the  same  form  as  in 
equation  (10),  namely 

A  =  -iPo^..3S(-iy(2/+l)Ah;(/:,,3rJj,(/:,3rJpXcos0j  for  0  <  r,  <  (13) 

1=0 

where  the  A  ^^e  coefficients  to  be  determined. 

The  radial  displacements  in  the  fluids  are  found  by  using 


to  yield 


M  =  — r~r'^P 
47tV^p 


1 


=  J^£(-l)'(2/+l)P,(cos0jh;(/:,,rJ  + A  ^h;(x) 

47rV'p,  S  I  dx  ax 
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„  = _ i ^ 

dr 


^(_i)'(2/+i)p^(cos9jh;(^^3r  )D;^j^(y) 
471  /  p3  ^  dy 


for  the  external  and  internal  fluids  respectively,  where  x  -  and  y  = 

To  find  the  displacements  and  stresses  in  the  shell  requires  the  solution  of  the  wave 
equation, 

clMv.u)-clWxVxu  =  ^.  (19) 

at 

The  displacement  is  derived  from  a  scalar  and  vector  potential. 


where 


m  =  -V4^  +  VxA 


cl, 


-Vx VxA  = 


1  a'A 
a  dl^  ■ 


Equation  (21)  is  of  the  same  form  as  equation  (6),  and  can  be  solved  as  in  Appendix  A 
to  give 


where 


^=vi/(r,e)T(0 


00 

V  =  S  P/  )[^/ j'  )  +  Fini(kL,2fr  )]• 


H;  is  the  spherical  Neumann  function  of  order  /.  Again  using  the  spherical  coordinate 
system,  separating  out  the  time  and  position  variables  and  noting  that  A  =  (o,0,  A^), 
equation  (22)  can  be  solved  (see  Appendix  B)  to  yield 

A,=a,(r,e,ct)).r(0  (25) 
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where 


1=0 

Returning  to  equation  (20),  the  displacements  in  the  shell  are 


ld\\l  1  d  (  . 


=  X  sin  0|-[£J,  (/cz„2?;  )+F'iT^i(ki,2rr  )]+-[GJ/  (^r';  )+'^/n,  (/c^r, )] 

1=0  I 


«$,2  =  0 


where  Tj  =  cos 0^,  z  =  and  T  =  kj-r^ . 
The  stresses  in  the  shell  are  (Faran,  1951) 


[r  r]  —  2p2Cy.2 


f7  ^  ,  ^‘^r.2 

- V.  M  H — r — 


r  ,  3Mb,  Mft,  1  9m_2 


r,^]  = 

^  rsin0  30  dr  r 


where  ct  and  g  are  Poisson’s  ratio  and  the  Lame  elastic  constant  of  the  shell  material, 
respectively.  Using 
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-  1^/2  ^  1  ^  (■  a  \  1 

.u  — \r  M  -  )H - - IsmS.Mfl,  )H - 

0r'‘  ’  rsin0  30^  ’  rsin0  9(j) 

=  2  P;  (cos  0,  )kl^  [eJ;  (k^^r, )  +  FiUi  (k^^r, )] 


these  become 


[rr]  =  2p242XP;(cos0,V^/:^%[£O;(^:i,2r,)+f;n,(/:,,2r,)]-;:22  E^^]i{z)+F,—n^{z) 


■  ^  2  ^  [Qj/  (^r^r  )  +  P^/ri/  (^T^r  )]  — ^  Q  “  ji  ^  ) 


=  +  EiYjiiz)  +  F,^n,iz) 

+  +  +  G,^j,(x)  +  //,^n,(x)  |  (38) 


(39) 

Substitution  of  equations  (5),  (10),  (12),  (13),  (16),  (18),  (28),  (30),  (37),  (38)  and  (39)  into 
equation  (4)  and  evaluation  at  r  =  a  and  r  =  b  yields  eight  simultaneous  equations 
from  which  the  coefficients  Q  needed  to  evaluate  the  scattered  pressure  can  be 
determined.  After  some  simplification  these  are  given  by 
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Yn 

Y12 

Y13 

Yi4 

^1 

0 

Y21 

Y  22 

Y  23 

Y24 

0 

Y26 

Y31 

Y  32 

Y  33 

Y34 

^3 

0 

Y41 

Y42 

Y43 

Y44 

0 

Y46 

Y  51 

Y52 

Y  53 

Y  54 

0 

0 

Y61 

Y62 

Y  63 

Y64 

0 

0 

Yn 

Yi2 

Yi3 

Yi4 

Yi5 

0 

Y21 

Y  22 

Y  23 

Y  24 

0 

Y26 

Y  31 

Y  32 

Y  33 

Y34 

Y  35 

0 

Y41 

Y42 

Y43 

Y44 

0 

Y46 

Y51 

Y52 

Y  53 

Y  54 

0 

0 

Y61 

Y62 

Y  63 

Y64 

0 

0 

The  coefficients  required  to  evaluate  equation  (40)  are  listed  in  Appendix  C.  The 
coefficients  D,  needed  to  calculate  the  pressure  in  the  filling  fluid  are  given  by 


Yn  Yi2  Yi3  Yi4  Yi5 

Y  21  Y  22  Y  23  Y  24  0  0 

Y  31  Y  32  Y  33  Y  34  Y  35  ^3 

Y41  Y42  Y43  Y44  0  0 

YSI  Y52  Y53  Y54  0  0 

kL,l^*(^L,l^p]  ^61  Y62  Y63  Y64  0  0 

kL3h;fkL3J'pl  Yi3  Yi4  Y  15  0 

^  '  ^Y21  Y22  Y23  Y24  0  Y  26 

Y31  Y32  Y33  Y34  Y  35  0 

Y41  Y42  Y43  Y  44  0  Y  46 

Y  51  Y  52  Y  53  Y  54  0  0 

Yei  Y62  Y63  Y64  0  0 


where  the  elements  of  the  determinants  are  again  as  listed  in  Appendix  C. 

Inserting  the  coefficients  of  equation  (40)  into  equation  (12)  yields  the  scattered 
pressure  at  the  listening  point  of  r  =  and  0  =  0^.  A  convenient  way  of  presenting 
this  pressure  is  to  relate  it  to  that  incident  upon  the  sphere,  in  the  form 


20 log  \rr 


If  0  =  180°  and  r,  is  sufficiently  large  so  that  spherical  spreading  applies,  i.e.  in  the  far- 
field,  then  this  term  equals  the  backscatter  target  strength,  defined  by 
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f 

r5  =  201og 


Ps 


le=180°,r=l 


Pi 


J 


(43) 


where  the  scattered  pressure  is  evaluated  1  m  from  the  centre  of  the  sphere  along  the 
axis  to  the  source.  In  the  comparisons  reported  below,  the  scattered  pressure  will  be 
given  as  the  target  strength  via  expression  (42)  with  6  =  180°  and  both  and  /; 
generally  set  to  1000  m  to  ensure  equivalence  with  equation  (43). 


2.2  Pulsed  Signals 

Until  now  the  point  source  has  been  considered  as  producing  continuous  sinusoidal 
waves.  However,  in  an  experiment  tone  bursts  or  chirped  signals  are  more  likely  to  be 
used.  The  scattered  pressure  for  these  two  cases  will  now  be  considered. 


2.2.1  Tone  Bursts 


Consider  first  a  tone  burst  of  central  frequency  /q  and  duration  At.  Defining 


.0  _27C/o 
'^L.X  ~ 


and 


rect(r)  = 


for  -y^<t<y2 

elsewhere 


we  have,  from  equation  (1), 


Pe 


-iklR 


R 


-  e‘^’''^‘’'rect 


—I  t-- 


At 


CJJ 


(44) 


(45) 


(46) 


which  has  the  Fourier  transform 


-Af/  p  p~^LX^ 

P.(f)^  f  -0 - 


sin[7t(/-/o)Af] 


(47) 

(48) 


as  the  frequency  components.  As  this  is  the  product  of  a  continuous  wave  mono¬ 
frequency  source  and  a  weighting  function,  the  scattered  pressure  frequency 
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components  may  be  represented  by  the  continuous  wave  solution  multiplied  by  the 
same  weighting  function.  Thus 


Psif)  =Psif) 


sin[7t(/-/o)Ai] 


(49) 


where  is  given  by  equation  (12).  Applying  the  inverse  Fourier  transform,  we  obtain, 
as  a  function  of  time. 


^.(0=1  Psif) 


4/-/o) 


(50) 


2.2.2  Chirps 


For  a  chirped  signal  of  start  frequency  /q,  stop  frequency  /i,  central  frequency  and 
duration  At,  equation  (46)  is  altered  to 


pM)  = 


P,e-''^- 


R 


K  i2n 

e 


4'+ 


2Al 


(If  R 
■'rect  —  t — 
\At\  c 


(51) 


where 


= 


(52) 


The  incident  signal  frequency  components  are 


Pff) 


P^e' 


% 


r/2 


2At 


dt. 


R 

The  integral  in  equation  (53)  is  evaluated  in  Appendix  E  to  give 

i27r[ 


2At 


dt 


At 


2(f-U) 


exp 


-i%At[l-ff 

(/i-/o) 


1+i 


erf 


(l-i). 


lAt 


f-f 


(/,-/o) 


(53) 


At 


t  +  fc-f 


Ah 


“/2 


Jj-y, 

(54) 
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The  scattered  pressure  frequency  components  are 


n(/) =?,(/)]_ 


^2 


i2d 


2A( 


dt 


(55) 


where  is  again  given  by  equation  (12).  As  a  function  of  time 


/>,(<)  =J_>,(/) 


i27C 


(/,-/)'+ 


(/|-/o)^2 
2  A/ 


dt 


>e 


(56) 


In  both  equations  (55)  and  (56)  the  integral  over  time  can  be  replaced  by  the  expression 
on  the  right  hand  side  of  equation  (54). 

In  both  equations  (50)  and  (56)  the  scattered  pressure  varies  in  amplitude  as  a 
function  of  time.  If  the  tone  burst  or  chirp  lasts  for  a  time  At,  the  backscatter  target 
strength  can  be  considered  as 


(57) 


(58) 


The  right  hand  side  of  equation  (58)  is  equivalent  to 


20  log 


1  rr+A; 


PM') 


2  ^ 
dt' 


(59) 


when  0  =  1 80°  and  r  and  are  large. 


It  follows  from  equation  (58)  that  the  target  strength  can  be  calculated  for  any  value 
of  t  corresponding  to  different  portions  of  the  return  signal.  The  most  intense 
reflections  will  be  the  specular  reflection  from  the  front  of  the  sphere,  and  the  first 
focused  return  from  the  rear  of  the  sphere.  For  6  =  1 80°  these  start  when 


and 


rp  +  r,.-  2b  At 


(60) 


11 


DSTO-RR-0020 


r^  +  r^-2b  2(b-a)  Aa  At 

t--P — +  + -  (61) 

Cj  C2  2 

respectively.  In  the  results  reported  below,  the  target  strength  of  the  sphere  will  be 
considered  as  the  naaximum  of  the  two  values  calculated  via  expression  (59)  at  these 
two  starting  times. 


2.3  Binary  Fluids 

In  practical  applications  of  the  above  theory  it  is  common  to  use  a  mixture  of  two 
liquids  for  the  interior  fluid.  For  a  binary  liquid  solution  the  longitudinal  sound  speed 
of  the  interior  fluid  becomes 


c.  - 


v^3  y 


where 


P3  = 


H.ap3,n  ^3.&P3.fc 


~  ^m,3fc^3,fc 


^3  “  ^3,a^i.a  ■*“  ^Xb^Xb' 


(62) 


(63) 

(64) 

(65) 


V3  ^  and  V3  ^  are  the  volumes  of  the  two  liquids,  P3  ^  and  P3  ^  are  their  densities,  3^ 

and  v^2b  molecular  velocities,  and  are  the  mole  fractions,  and 

and  M^h  ^re  the  molecular  masses.  In  general,  the  densities  and  velocities  are  a 
function  of  temperature,  so  the  speed  of  sound  of  the  interior  fluid  will  likewise  be 
temperature  dependent. 


3.  Physical  Data 

Consider  a  spherical  shell  filled  with  a  fluid  or  mixture  of  two  fluids  and  placed  in  the 
sea.  If  a  small  sound  projector  directs  a  continuous  sinusoidal,  tone  burst  or  chirped 
signal  at  it,  the  reflected  signal  picked  up  by  a  hydrophone  in  the  same  or  another 
location  will  depend  on  a  large  number  of  factors.  Parameters  from  the  incident  sound 
include  the  distance  from  the  source  to  the  centre  of  the  sphere,  the  frequency  for  a 
continuous  signal,  the  frequency  and  pulse  length  for  a  tone  burst,  and  the  start  and 
stop  frequencies  and  pulse  length  for  a  chirped  signal.  The  distance  from  the  centre  of 
the  sphere  and  angular  position  of  the  receiver  are  other  factors,  as  are  a  number 
associated  with  the  reflecting  sphere.  These  include  the  diameter,  thickness,  and 
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material  out  of  which  the  shell  is  made.  The  latter  factor  determines  the  density  of  the 
shell,  as  well  as  Poisson's  ratio  for  the  material,  and  the  velocity  of  longitudinal  and 
shear  waves  in  it.  The  temperature,  salinity  and  depth  of  the  propagation  path  set  the 
density  and  sound  velocity  of  the  fluid  surrounding  the  sphere.  The  temperature  and 
the  filling  fluid  or  mixture  of  fluids  chosen  affect  the  interior  fluid's  density  and  sound 
velocity. 

Clearly,  when  long  projector-to-sphere  ranges  are  used,  sound  absorption  in  the 
water  may  be  significant.  As  this  will  not  affect  a  comparison  of  different 
sphere /filling  fluid/signal  type  parameters,  it  has  not  been  included  in  the  calculations 
reported.  Nor  has  absorption  within  the  sphere  or  its  filling  fluid,  as  the  ranges 
considered  are  so  small. 

Table  1  lists  the  parameters  for  the  shell  materials  studied.  For  the  fluids  the  density 
and  sound  velocity  depend  on  a  number  of  factors,  as  noted  above,  so  equations  are 
needed  to  express  these  dependences.  Millero  and  Poisson  (1981)  give  the  density  of 
sea  water  (in  kgm“^)  as 


(  \ 


(66) 


where 

p°  =  999.842594  +  6.793952  x  lO’^f  -  9.095290  x  + 1.001685  x  10^ 
-1.120083  X 10"^  +  6.536336  x  10"^  f ^  +  (8.24493  x  10'^  -  4.0899  x  10“^  t 

+7.6438  X 10"^  -  8.2467  X 10"^  +  5.3875  X 10"^  )S  +  (-5.72466  X 10-3 
+1.0227  X 10-^  f- 1.6546  X 10"^  f  2  +4.8314x10-^52  (67) 


K  =  19652.21  +  148.4206f  -  2.327105t2  + 1.360477  x  IQ-^  -  5.155288  x  IQ-^ 

+S(54.6746  -  0.603459f  + 1.09987  x  10-^  -  6.1670  x  10-^  f ^  ^  +  (7.944  x  lO'^ 

+1.6483  X 10-2  f-  5.3009  x  10-^f2)+  p[3.239908  + 1.43713  x  10-3  f  +  j  15032  x  10-^  f  2 
-5.77905  X 10-2 P  +  s(2.2838  x  10-3  -  ^  0931  x  lO'^ t  - 1.6078  x  10"^ ) 

’2 [8.50935  X 10-3  -  6.12293  x  10-^ t  +  5.2787  x  10-^ f 2 


+5/2(1.91075x10-^) 

+s(-9.9348  X 10-2  ^  2.0816  x  10-^  f  +  9.1697  x  10-^°  ) 


(68) 


and 
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P^\p^gdy.  (69) 

P  is  the  gauge  pressure  in  bars,  with  P  =  0  at  the  water  surface  under  1  atm  of  air 
pressure.  A  rough  rule  of  thumb  is  that  the  pressure  increases  by  1  bar  for  every  10  m 
of  water  depth.  S  is  the  salinity  in  %o,  t  is  the  temperature  in  °C,  d  is  the  depth  in  m 
and  g  is  the  acceleration  due  to  gravity.  The  speed  of  sound  (in  ms  ')  is  given  as 

c,  =1412  +  3.2k  +  1.19S  +  0.0167fl!  (70) 

by  Francois  and  Garrison  (1982). 


Table  1: 


Material 

Density 
(kg  m'^) 

Poisson's  ratio 

Velocity  of  longitudinal 
waves  (ms ') 

Velocity  of  shear  waves 
(ms^) 

Aluminium 

2700 

0.355 

6420 

3040 

Brass 

8600 

0.374 

4700 

2110 

Nickel 

8900 

0.336 

6040 

3000 

Pyrex  glass 

2320 

0.17 

5640 

3280 

Steel 

7910 

0.29 

5790 

3100 

The  above  data  were  obtained  from  Gray  (1972). 


To  achieve  a  high  target  strength  in  the  backscatter  direction,  a  filling  fluid  with  a 
refractive  index  between  1.5  and  2.0  is  needed  to  obtain  the  necessary  thin  lens 
focusing.  This  is  most  easily  accomplished  by  mixing  two  liquids  of  different  sound 
speeds,  one  above  and  the  other  below  the  required  speed.  Other  considerations  are 
cost  and  safety  of  the  liquids.  Some  of  the  potential  liquids,  such  as  CCI4,  are 
hazardous.  Initially  Freon-llS’^'^  and  ethanol,  mixed  in  varying  ratios,  were  studied  as 
the  filling  liquid,  as  they  are  cheap,  relatively  safe  and  readily  available.  However  due 
to  expected  difficulties  in  obtaining  Freon-1  Fluorinert  FC-72'^'^  and  n-hexane 
have  also  been  studied.  From  Gray  (1972)  the  density  and  speed  of  sound  in  ethanol 
are 


P3^  =785.06-0.8591(f-25)-5.6xl0^(f-25)^  -  5xl0-®(r-25)^  for  10<t<40 

(71) 


and 

=  1232.1 3- 3.5739/ -b 2.39 xlO'^t^  -b 2.65 x  10'^^  - 2.4 x  10"^/'’  for  0</<58. 

(72) 

Fitting  curves  to  the  data  of  Leader  et  al.  (1969),  the  density  and  speed  of  sound  in 
Freon-11 3^”^  are 
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p3^  =  1614.5-2.04?  for0<?<40 


(73) 


and 


j786.l-3.345?  for  0<?<20 

[778.7-2.975?  for  20<?<40' 


(74) 


Fitting  a  curve  to  the  data  of  Timmermans  (1950),  the  density  of  n-hexane  is 

p3„  =  676.891-0.846?-l.lxl0■^?^  (75) 

Kaye  and  Laby  (1973)  quote  the  speed  of  sound  in  it  at  30°C  as  1060  ms"‘,  so  using 
the  relation 


(76) 


where  B  is  the  bulk  modulus,  the  speed  of  sound  in  n-hexane  can  be  given  as 


650.5x1060' 


P3.a 


(77) 


3M's  Fluorinert™  Liquids,  Product  Manual  gives  the  density  of  Fluorinert  FC-72™  as 

P3^  =1740-2.61?.  (78) 

From  the  data  for  the  speed  of  sound, 

C3^  =582-2.8?.  (79) 

In  the  results  reported  below,  comparison  will  be  made  with  an  air-filled  sphere.  At 
sea  level,  Kaye  and  Laby  (1973)  give  the  density  of  air  for  50%  relative  humidity  as 

^  1.29307xl01.3-0.2443p„  ^  1. 000028 x  (1  +  0.00367?)  (80) 

where  the  saturated  vapour  pressure  p„can  be  found  by  fitting  a  curve  to  their  data  to 
yield 

p,,  =0.6116  +  4.365xl0"'?  +  1.567xl0"'?' +1.74x10"'?' +5.5x10"' ?\  (81) 
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The  speed  of  sound  in  air  for  sound  waves  of  frequency  /  is  (Kneser,  1931) 


(82) 


where  the  zero-frequency  speed  can  be  found  by  substituting  into  the  equation  of 
Miller  (1937)  to  yield 


c.  =  331.46xJ  !  +  ■ 


273 


X 


,  0.5p„  ri.324  18.01534 

1 - ^x 


101.3 


U.403  28.9644 


(83) 


The  infinite-frequency  speed  c^exceeds  this  by  0.12  ms  '  at  20°C  (Harris,  1971).  The 
relaxation  frequency  at  50%  relative  humidity  is  (Piercy,  1969) 


fr 


3.05x10' 


50p, 


.1,3 


0.5p,,4-101.325 


(84) 


4.  Computational  Considerations 

The  summation  in  equation  (12)  is  theoretically  over  an  infinite  number  of  terms, 
although  in  practice  the  summation  was  halted  when  the  relative  difference  of  the 
scattered  pressure  after  summation  to  the  /th  and  l-lth  terms  was  less  than  0.01%  and 
after  summation  to  the  /-1th  and  /-2th  terms  was  also  less  than  0.01%.  To  ensure 
convergence  care  had  to  be  taken  in  calculating  the  spherical  Bessel  functions  and  their 
derivatives.  Details  of  the  techniques  used  are  given  in  Appendix  D. 

As  a  check  on  the  formulae,  computational  techniques  and  the  computer  programs 
used  to  calculate  the  results  presented  in  the  next  section,  the  magnitude  of  the  form 
function,  defined  by 


2 


/=o 


(85) 


was  calculated  for  an  iron  shell  of  20  cm  diameter,  20  mm  wall  thickness,  longitudinal 
and  shear  wave  sound  speeds  of  5950  and  3240  ms”',  respectively,  and  Poisson  ratio 
of  0.292.  This  shell  was  immersed  and  filled  with  water  of  density  1000  kgm”^  and 
sound  speed  1410  ms”’.  The  results  are  displayed  in  Fig.  2  along  with  those  of 
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Hickling  (1964)  for  a  similar  set  of  parameters.  The  agreement  is  generally  good, 
although  there  are  some  variations  attributed  to  two  sources.  The  magnitude  and 
position  of  some  of  the  peaks  and  troughs  do  not  align  perfectly,  due  to  a  small 
inconsistency  in  Hickling's  shear  and  Lame  elastic  constants  and  longitudinal  and 
shear  sound  speeds  for  iron.  In  addition,  the  narrow  spikes  are  sensitive  to  the 
evaluation  frequency,  which  was  not  identical  for  both  data  sets. 


5.  Results 

Fig.  3  shows  the  calculated  absolute  value  of  the  scattered  pressure  field  in  a  4  m  by  4 
m  section  centred  about  a  stainless  steel  spherical  shell  of  20  cm  diameter  and  0.8  mm 
wall  thickness.  The  sphere  is  indicated  by  the  white  circle  and  is  filled  with  a  mixture 
of  34%  Fluorinert  FC-72’’‘''^  and  66%  n-hexane  by  volume  and  deployed  at  shallow 
depths  in  sea  water  of  20°C  and  35  %o  salinity.  The  ratio  of  filling  fluids  was  chosen  to 
give  a  refractive  index  of  1.8  at  20°C.  The  projector  is  located  1000  m  to  the  left  of  the 
origin  and  1  m  from  it  (the  projector)  the  amplitude  was  unity.  The  scattered  pressure 
amplitudes  are  relative  to  this  value.  The  projector  emits  continuous  waves  of  100 
kHz  frequency.  The  field  shows  strong  focusing  in  the  forward  scattering  direction,  as 
well  as  the  presence  of  weak  intensity  regions  and  nodes.  There  is  also  some  focusing 
in  the  backscatter  direction. 

Fig.  4  shows  the  scattered  pressure  amplitude  field  with  the  same  sphere  and 
environment  parameters.  The  projector  emits  continuous  waves  of  200  kHz.  Note  the 
narrower  focusing  in  the  backscatter  direction. 

Fig.  5  displays  the  total  (incident  plus  scattered)  pressure  amplitude  field  in  a  section 
through  the  centre  of  the  sphere,  in  the  vicinity  of  the  sphere,  whose  inner  and  outer 
walls  are  shown  by  the  black  circles.  The  parameters  are  as  for  Fig.  3.  The  waves  to 
the  left  show  what  may  be  spherical  spreading  but  with  a  radius  of  curvature  slightly 
larger  than  that  of  the  sphere.  To  the  right  of  the  sphere  the  nodes  and  antinodes  form 
a  pattern  radial  from  a  point  somewhat  to  the  left  of  the  centre  of  the  sphere.  The  wave 
pattern  inside  the  sphere  is  complex.  A  set  of  nodes  and  antinodes  occurs  both  radially 
and  circumferentially,  with  the  greatest  intensity  along  the  axis  from  the  projector. 

Fig.  6  shows  the  total  pressure  amplitude  field  under  the  same  conditions  for 
continuous  waves  of  200  kHz.  The  features  are  similar  to  that  for  Fig.  5,  but  on  a 
smaller  scale  commensurate  with  the  higher  frequency. 

Fig.  7  overlays  four  sets  of  calculated  backscattered  target  strength  values  for  a 
projector  and  receiver  collocated  1000  m  from  a  sphere.  The  sphere  is  of  stainless  steel 
of  20  cm  diameter  and  0.8  mm  wall  thickness  and  is  filled  with  a  mixture  of  68%  Freon- 
113^“  and  32%  ethanol  by  volume  and  deployed  at  shallow  depths  in  sea  water  of  20°C 
and  35  %o  salinity.  The  ratio  of  filling  fluids  was  chosen  to  give  a  refractive  index  of  1.8 
at  20°C. 
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The  wildly  varying  line  is  the  target  strength  as  a  function  of  frequency  for  steady 
state  conditions.  The  target  strength  was  calculated  in  100  Hz  steps  from  100  Hz  to 
500  kHz.  For  very  low  frequencies  the  sphere  is  small  compared  to  the  wavelength  of 
sound  in  water,  and  very  little  sound  is  reflected.  With  increasing  frequency  the  target 
strength  increases,  but  fluctuates  wildly  with  small  changes  in  the  frequency.  Beyond 
450  kHz  it  takes  on  an  oscillatory  nature. 

The  target  strength  as  a  function  of  frequency  in  1  kHz  steps  to  450  kHz  is  shown  for 
both  100  and  500  ps  tone  bursts.  The  result  for  500  ps  chirps  of  10  kHz  bandwidth 
centred  on  frequencies  taken  in  1  kHz  steps  to  450  kHz  is  also  shown.  Although  there 
are  small  variations  between  the  last  three  curves,  especially  below  50  kHz  and  above 
400  kHz,  these  are  of  a  minor  nature. 

Fig.  8  shows  the  fine  detail  of  the  steady  state  data  of  Fig.  7  between  100  and  120  kHz, 
calculated  in  25  Hz  steps.  It  is  compared  with  data  in  100  Hz  steps.  The  overall  result 
is  the  same,  although  several  narrow  spikes  are  evident  in  the  former  set  of  data.  The 
curves  indicate  that  most  of  the  major  target  strength  variations  are  somewhat 
oscillatory,  with  a  period  of  approximately  1.6  kHz.  Overlaid  on  this  is  a  fine  structure 
with  a  width  of  less  than  25  Hz. 

The  target  strength  as  a  function  of  frequency  in  1  kHz  steps  from  100  to  120  kHz  is 
also  shown  in  Fig.  8  for  100  ps  tone  bursts.  One  curve  shows  the  results  based  on 
steady  state  data  calculated  in  25  Hz  steps;  the  other  shows  the  result  when  the  steady 
state  data  has  been  calculated  in  100  Hz  steps.  The  results  are  almost  identical, 
indicating  that  the  fine  structure  has  little  effect  for  100  ps  tone  bursts. 

The  complexity  of  the  target  strength  as  a  function  of  frequency  for  both  continuous 
and  pulsed  waves  is  due  to  the  existence  of  resonances  of  several  types,  which  occur  in 
the  fluid  and  the  shell.  At  low  frequencies  "Franz"  or  "creeping"  waves  circumnavigate 
the  shell  in  the  external  fluid.  The  eigenfrequencies  depend  upon  the  shape  and  size  of 
the  shell  and  correspond  to  constructive  interference  between  the  specularly  reflected 
waves  and  those  circumnavigating  the  sphere.  "Rayleigh"  and  "whispering  gallery" 
waves,  which  depend  on  the  composition  of  an  elastic  target,  occur  at  frequencies  for 
which  multiple  half-integral  wavelengths  correspond  to  the  circumference,  thereby 
producing  standing  waves  on  the  surface  which  radiate  back  into  the  fluid.  These  are 
applicable  to  solid  elastic  targets.  "Lamb"  waves  are  the  analogue  for  the  elastic  shell, 
with  extensional  and  flexural  resonances  occurring.  They  depend  on  the  composition 
of  the  shell.  Other  resonances  occur  when  the  multiple  half-integral  wavelength 
associated  with  the  compressional  waves  in  the  shell  equals  the  shell  thickness. 

Figs  9  and  10  display  the  target  strength  as  a  function  of  frequency  to  450  kHz  for 
refractive  indices  of  1.5, 1.6, 1.7, 1.8  and  1.9.  In  each  diagram  the  projector  and  sphere 
are  deployed  at  shallow  depths  in  sea  water  of  20°C  and  35 %o  salinity.  The  projector 
and  receiver  are  located  1000  m  away  and  the  former  emits  100  ps  tone  bursts.  The 
spherical  shell  is  made  of  stainless  steel  and  is  20  cm  in  diameter  and  of  0.8  mm  wall 
thickness.  In  Fig.  9  the  shell  was  filled  with  a  mixture  of  Freon-llS"^*^  and  ethanol;  in 
Fig.  10  the  filling  liquids  were  Fluorinert  FC-72’'’^  and  n-hexane.  The  ratios  of  the  two 
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sets  of  liquids  were  chosen  to  give  the  same  velocities  of  sound,  but  the  densities  were 
not  identical.  In  general,  the  target  strengths  with  the  Fluorinert  FC-72™  and  n-hexane 
combination  are  less  than  for  the  Freon-113™  and  ethanol  combination,  but  the  overall 
shape  of  each  pair  of  curves  is  similar. 

Figs  11  and  12  show  the  effect  of  shell  material  composition  on  the  target  strength  for 
100  [Ls  tone  bursts  of  frequencies  to  450  kHz.  Again  the  projector  and  receiver  were 
located  1000  m  from  a  20  cm  diameter  shell  of  0.8  mm  wall  thickness  in  shallow  depths 
of  sea  water  of  20°C  and  35  %o  salinity.  In  Fig.  11  the  shell  was  filled  with  68%  Freon- 
113'^“  and  32%  ethanol  by  volume;  in  Fig.  12  the  liquids  were  34%  Fluorinert  FC-72^^ 
and  66%  n-hexane  by  volume.  In  each  case  the  refractive  index  was  1.8.  The  structure 
of  the  curves  is  similar  for  each  pair  of  liquids,  but  again  the  target  strengths  with  the 
Fluorinert  FC-72™  and  n-hexane  combination  are  less  than  for  the  Freon-113™  and 
ethanol  combination.  Generally,  glass  and  aluminium  yield  substantially  larger  target 
strengths  than  for  the  other  materials,  although  stainless  steel  and  nickel  are  more 
likely  to  be  used  as  targets  in  practice  because  of  their  greater  durability. 

Concentrating  on  stainless  steel  as  the  shell  material.  Fig.  13  shows  the  effect  of 
varying  the  shell  diameter.  All  other  parameters  are  as  for  Fig.  12.  As  expected,  the 
target  strength  increases  with  diameter.  The  oscillatory  behaviour  evident  in  Fig.  7  is 
very  pronounced  for  the  10  cm  diameter  sphere. 

Keeping  to  a  stainless  steel  shell.  Figs  14  and  15  display  the  effect  of  varying  the  wall 
thickness.  All  other  parameters  are  as  for  Figs  11  and  12,  respectively.  Again,  the 
Fluorinert  FC-72™  -  n-hexane  combination  results  in  lower  target  strengths  than  the 
Freon-113™  -  ethanol  combination.  For  very  thin  walls,  which  are  not  physically 
practical,  the  target  strength  increases  relatively  smoothly  with  frequency,  and  is  less 
than  that  for  thicker  walls  at  low  frequencies.  As  the  thickness  increases,  the  target 
strength  rises  to  a  maximum  with  frequency,  then  declines.  The  frequency  and  height 
of  the  maximum  decreases  with  increasing  wall  thickness.  As  the  wall  thickness 
increases  the  oscillatory  behaviour  referred  to  above  starts  to  appear,  and  begins  at 
progressively  lower  frequencies.  The  curve  for  1.6  mm  wall  thickness  indicates  that 
this  oscillatory  response  does  not  continue  for  all  higher  frequencies,  but  only  occurs 
over  a  limited  frequency  range,  before  dying  out  and  returning  to  a  more  regular 
response. 

Figs  16  and  17  show  how  the  refractive  index  varies  for  changes  in  the  percentage  by 
volume  of  Freon-113™  and  ethanol,  and  Fluorinert  FC-72™  and  n-hexane,  respectively. 
The  curves  are  plotted  for  waters  of  35  %o  salinity  and  water  (and  filling  liquid) 
temperatures  of  0,  5, 10, 15,  20,  25,  30  and  35°C.  Both  sets  of  curves  show  sensitivity  of 
the  refractive  index  to  temperature,  and  since  Figs  9  and  10  in  turn  indicate  sensitivity 
of  the  target  strength  to  refractive  index,  it  can  be  expected  that  with  all  other 
parameters  constant,  the  target  strength  will  change  with  temperature. 
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Freon-113™  and  32%  ethanol  by  volume  and  deployed  at  shallow  depths  in  sea  water  of20°C  and  35%o  salinity.  The  target  strength  is  shown 
for  the  shell  materials  stainless  steel,  aluminium,  nickel,  brass  and  pyrex  glass  when  the  incoming  waves  are  100  |is  tone  bursts. 
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Fluorinert  FC-72™  and  66%  n-hexane  hy  volume  and  deployed  at  shallow  depths  in  sea  water  of  20°C  and  35%o  salinity.  The  target  strength 
shown  for  shell  diameters  of  10,  20  and  30  cm  when  the  incoming  waves  are  100  jxs  tone  bursts. 
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Freon-113™  and  32%  ethanol  by  volume  and  deployed  at  shallow  depths  in  sea  water  of  20°C  and  35%o  salinity.  The  target  strength  is  shown 
for  wall  thicknesses  of  0.01,  0.2,  0.4,  0.6, 0.8, 1.0, 1.2, 1.4  and  1.6  mm  when  the  incoming  waves  are  100  |is  tone  bursts. 
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ire  1 7:  Refractive  index  as  a  function  of  the  percentage  of  Fluorinert  FC-72™  in  a  Fluorinert  FC-72™  -  n-hexane  mixture,  wi 
exterior  water  of  35%o  salinity  are  maintained  at  temperatures  of 0,  5, 10, 15,  20,  25,  30  and  35°C. 
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Figs  18  and  19  show  how  the  target  strength  alters  for  100  )is  tone  bursts  of 
frequencies  up  to  450  kHz,  for  temperatures  of  10,  20  and  30°C.  Again  the  projector 
and  receiver  were  located  1000  m  from  a  20  cm  diameter  stainless  steel  shell  of  0.8  mm 
wall  thickness  in  shallow  depths  of  sea  water  of  35%o  salinity.  In  Fig.  18  the  shell  was 
filled  with  68%  Freon-113^'^  and  32%  ethanol  by  volume;  in  Fig.  19  the  liquids  were 
34%  Fluorinert  FC-72''''^  and  66%  n-hexane  by  volume.  Below  50  kHz  temperature 
changes  produce  only  a  small  change  in  the  target  strength,  but  for  higher  frequencies 
the  changes  get  progressively  larger,  especially  in  going  from  10  to  20®.  Of  course 
smaller  temperature  variations  will  alter  the  target  strength  by  correspondingly  less. 

Fig.  20  compares  the  target  strengths  for  100  ps  tone  bursts  of  frequencies  to  450  kHz, 
for  a  20  cm  diameter  stainless  steel  spherical  shell  of  0.8  mm  wall  thickness  deployed  at 
shallow  depths  in  sea  water  of  20°C  and  35 %o  salinity.  The  filling  fluids  were  68% 
Freon-113'^^  and  32%  ethanol  by  volume,  34%  Fluorinert  FC-72'^'^  and  66%  n-hexane  by 
volume,  air  of  50%  relative  humidity,  and  water  of  35 %o  salinity.  The  target  strength 
for  the  air-filled  sphere  is  fairly  constant  and  in  close  agreement  to  the  formula  for 
specular  reflection  from  a  sphere  whose  size  is  large  compared  to  the  wavelength,  viz. 

I  2 

75  =  10  logy.  (86) 


Above  100  kHz  the  target  strength  of  the  water-filled  sphere  shows  periodic  changes  of 
about  10  dB,  with  a  period  in  the  frequency  domain  of  approximately  90  kHz.  In 
general  the  values  are  higher  than  for  the  air-filled  sphere.  However,  the  Freon-113™  - 
ethanol  and  Fluorinert  FC-72^“  -  n-hexane  filled  spheres  have  even  larger  target 
strengths,  which  vary  relatively  smoothly  with  frequency.  In  comparison  to  the  air- 
filled  sphere,  there  is  a  gain  in  excess  of  10  dB  in  target  strength  at  all  frequencies 
studied.  As  noted  previously,  the  performance  of  the  two  pairs  of  liquid  combinations 
is  similar,  but  not  identical. 


Although  the  theory  developed  is  not  restricted  to  the  far  field,  the  results  of  Fig.  7  to 

20  were  based  on  the  projector  and  receiver  being  1000  m  from  the  centre  of  the  sphere. 
Keeping  the  projector  1000  m  from  the  centre  of  the  sphere,  but  varying  the  distance  of 
the  receiver  out  along  the  axis  from  the  centre  of  the  sphere  towards  the  projector.  Fig. 

21  shows  the  variation  in  apparent  target  strength  for  a  20  cm  diameter  stainless  steel 
shell  of  0.8  mm  wall  thickness  deployed  at  shallow  depths  in  sea  water  of  20°C  and 
35 %o  salinity  when  filled  with  34%  Fluorinert  FC-72^^  and  66%  n-hexane  by  volume. 
At  100,  200,  300  and  400  kHz  the  farfield  condition 


i2by 

^L.l 


(87) 


is  satisfied  at  approximately  2.7,  5.3,  8  and  10.7  m,  respectively.  Examination  of  the 
curves  in  Fig.  21  shows  constancy  of  the  apparent  target  strength  occurs  at  ranges 
slightly  beyond  these  at  the  respective  frequencies.  Below  approximately  200  kHz,  the 


51 


DSTO-RR-0020 


apparent  target  strength  increases  with  range;  at  higher  frequencies  it  generally 
decreases  with  range. 

Figs  7  to  21  have  shown  intensities  in  the  frequency  domain.  Figs  22  to  25  show  the 
reflected  signal  in  the  time  domain,  where  0  ms  corresponds  to  the  arrival  of  the  first 
portion  of  the  reflected  signal  at  the  receiver.  The  pressure  is  relative  to  a  value  of  1  at 
the  projector.  In  each  diagram  the  waves  impinge  upon  a  20  cm  diameter  stainless 
steel  shell  of  0.8  mm  wall  thickness  deployed  at  shallow  depths  in  sea  water  of  20°C 
and  35  %o  salinity  filled  with  68%  Freon-113™  and  32%  ethanol  by  volume.  In  Fig.  22 
the  signal  emitted  by  the  projector  was  a  200  kHz  tone  burst  of  100  |xs  duration.  The 
most  intense  signal  corresponds  to  the  focused  reflection  from  the  rear  of  the  sphere, 
which  arrives  0.44  ms  after  the  front  surface  specular  reflection.  Subsequent  reflections 
are  related  to  multiple  internal  reflections  within  the  sphere.  Fig.  23  shows  the 
equivalent  signal  for  a  tone  burst  lasting  500  ps.  Due  to  the  long  duration  of  the  tone 
burst  there  is  overlap  of  signals  returning  from  different  parts  of  the  sphere,  but  that 
from  the  rear  surface  still  dominates. 

In  Figs  24  and  25  a  chirped  signal  of  10  kHz  bandwidth  centred  on  200  kHz  was 
used.  In  Fig.  24  the  duration  of  the  chirp  was  100  ps  and  although  the  signal  from  the 
rear  surface  dominates,  a  later  multiply-reflected  signal  is  almost  as  large.  Fig.  25 
shows  that  a  500  ps  chirp  results  in  much  interference  of  returns  from  different  parts  of 
the  sphere.  The  first  rear  surface  reflection  and  a  somewhat  later  multiply-reflected 
signal  are  the  most  prominent  returns. 


6.  Conclusions 

Selection  of  the  appropriate  materials  to  make  a  fluid-filled  spherical  shell  target  of  a 
required  target  strength  is  no  easy  matter,  for  not  only  is  the  target  strength  sensitive  to 
a  number  of  material  parameters,  but  with  a  given  selection  of  materials,  the  target 
strength  is  sensitive  to  both  the  form  of  the  incoming  sound  waves  and  the 
environment. 

Consider  first  the  incoming  sound  waves.  Fig.  7  shows  that  the  target  strength  varies 
greatly  for  small  changes  in  the  frequency  if  continuous  waves  are  used. 
Experimentally  the  implication  is  that  drift  in  the  frequency  will  lead  to  a  large  change 
in  the  target  strength.  Fortunately,  these  variations  can  be  largely  eliminated  by  the 
use  of  pulsed  signals,  which  effectively  average  over  many  frequencies.  However, 
there  are  some  target  strength  variations  with  pulse  length  below  about  50  kHz,  or  in 
the  region  of  oscillatory  target  strength  behaviour. 
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with  68%  Freon-113™  and  32%  ethanol  by  volume  and  deployed  at  shallow  depths  in  sea  water 
incoming  waves  are  100  jis  tone  bursts. 
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sphere  to  the  projector  and  is  0.2,  0.4,  0.6,  0.8, 1.0, 1.2, 2, 3, 5, 10  and  1000  m  from  the  centre  of  the  sphere. 
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Figure  22:  The  reflected  signal  arriving  at  a  receiver  co-located  with  a  projector  1000  mfrom  a  stainless  steel  spherical  shell  of  20  cm  diameter 
and  0.8  mm  wall  thickness  filled  with  68%  Freon-113'^**  and  32%  ethanol  by  volume  and  deployed  at  shallow  depths  in  sea  water  of20°C  and 
35%o  salinity.  The  projector  emitted  a  100  \is  tone  burst  of 200  kHz,  and  0  ms  corresponds  to  the  arrival  of  the  first  part  of  the  return  signal. 
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Figure  23:  The  reflected  signal  arriving  at  a  receiver  co-located  with  a  projector  1000  m  from  a  stainless  steel  spherical  shell  of  20  cm  diameter 
and  0.8  mm  wall  thickness  filled  with  68%  Freon-113'^**  and  32%  ethanol  by  volume  and  deployed  at  shallow  depths  in  sea  water  of  20°C  and 
357oo  salinity.  The  projector  emitted  a  500  |xs  tone  burst  of 200  kHz,  and  0  ms  corresponds  to  the  arrival  of  the  first  part  of  the  return  signal. 
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Figure  25:  The  reflected  signal  arriving  at  a  receiver  co-located  with  a  projector  1000  m  from  a  stainless  steel  spherical  shell  of  20  cm  diameter 
and  0.8  mm  wall  thickness  filled  with  68%  Freon-113^*'^  and  32%  ethanol  by  volume  and  deployed  at  shallow  depths  in  sea  water  of  20°C  and 
35%o  salinity.  The  projector  emitted  a  10  kHz  bandwidth  chirp  centred  on  200  kHz  and  of 500  |is  duration,  and  0  ms  corresponds  to  the  arrival  of 
the  first  part  of  the  return  signal. 
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Now  consider  the  material  parameters  which  will  lead  to  a  large  target  strength 
above  100  kHz.  A  shell  of  aluminium  or  pyrex  glass  would  optimise  the  target 
strength,  but  as  these  are  soft  or  brittle  materials,  stainless  steel  might  be  preferable  as 
it  is  more  robust.  As  might  be  expected,  increasing  the  diameter  of  the  sphere  will  give 
a  larger  target  strength.  A  thin  shell  wall  will  also  lead  to  a  large  target  strength,  but 
may  be  too  difficult  to  manufacture  and  will  be  easily  damaged.  A  thickness  of 
0.8  mm  might  be  classed  as  a  reasonable  compromise.  Figs  14  and  15  indicate  that  the 
target  strength  is  quite  sensitive  to  the  thickness,  so  variations  in  the  wall  thickness 
need  to  be  avoided.  In  particular,  a  sphere  made  by  welding  together  two  spun 
hemispheres  is  not  appropriate,  as  the  process  of  spinning  leads  to  a  variable  wall 
thickness. 

Choice  of  the  optimum  combination  of  fluids  with  which  to  fill  the  shell  depends  on 
the  temperature.  Variations  of  a  few  degrees  will  make  only  a  modest  change  to  the 
target  strength,  but  a  combination  chosen  for  warm  waters,  say  off  northern  Australia, 
will  behave  quite  differently  in  cool  waters  off  southern  Australia.  For  frequencies 
below  200  kHz  a  mixture  of  Freon-113™  and  ethanol,  or  Fluorinert  FC-72™  and  n- 
hexane  will  work  quite  well  if  the  ratio  of  liquids  is  chosen  to  give  a  refractive  index 
between  1.8  and  1.9.  At  frequencies  above  250  kHz  the  same  liquids  mixed  to  give  a 
refractive  index  between  1.6  and  1.7  will  produce  high  target  strengths. 

In  general,  a  well  selected  fluid-filled  sphere  can  give  a  large  target  strength,  but  it 
will  be  somewhat  sensitive  to  temperature  and  hence  refractive  index  of  the  filling 
fluids.  Consequently,  the  water  temperature  at  the  site  at  which  it  is  to  be  used  should 
be  measured  so  the  target  strength  can  be  calculated.  Better  still,  the  target  strength 
should  actually  be  measured  at  the  site. 

If  it  is  difficult  to  manufacture  a  shell  with  a  uniform  wall  thickness,  consideration 
could  be  given  to  spherical  shells  composed  of  rubber-coated  metals  or  non-metallic 
materials,  such  as  plastics.  It  is  possible  that  the  target  strengths  may  be  less  sensitive 
to  variations  in  their  shell  parameters.  Because  of  their  significant  viscosity,  the  above 
theory  is  not  valid,  so  they  were  not  compared  in  this  study.  The  required  theory  is  an 
extension  of  that  presented  here. 
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Appendix  A 

To  solve  for  P  in  equation  (6)  the  spherical  coordinates 

;c=rsin0cos(l)  (Al) 

=  rsinSsincj)  (A2) 

z  =  rcos0  (A3) 

are  used  and  it  is  noted  that  P  has  no  (j)  dependence.  Hence  equation  (16)  can  be 
rewritten  as 


AAf  _J_Afsin0— 1  =  — — 

dr  j  sin  0  30 1  307  c^, 


Solution  of  equation  (A4)  proceeds  by  separating  out  the  time  and  position 
dependence  of  P  as 

P  =  p{r,Q).T{t)  (A5) 

and  substituting  this  into  equation  (A4).  Hence 

I  df  2dp^  1  37  .  ^dp^  1  d^T 

r^pdr\  dr)  sin 0  30  V  30  j  c^jT  dt 

The  left  and  right  hand  sides  of  equation  (A6)  are  independent  and  can  now  be  set 
equal  to  a  constant  -7,,/  so  the  right  hand  side  becomes 


for  which  the  solution  is 


■+/cLci7  =  o 


From  the  left  hand  side  of  equation  (A6), 


i  A f V  . J„..  A 1 + /tj, p = 0 

r^dr\  dr)  r‘sin0  30v  dO) 


which  can  be  solved  by  a  separation  of  variables.  Substituting 

p(r.0)  =  /?(/-).  7e) 


(AlO) 
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into  equation  (A9)  yields 

\  d  (  2dR\  ,2  2  1  d  (  .  ^dY\ 

R  dr  \  dr )  y  sin  6  (i6  V  dQ  J 


(All) 


Again  the  left  and  right  hand  sides  of  equation  (All)  are  independent  and  are  set  equal 
to  /(/  + 1).  The  right  hand  side  becomes 


1 


sin0 


.dY  .  .d^Y 

cos  9 - i-sm6 — ^ 

de  dQ^ 


+  1{1+1)Y  =  0. 


which  can  be  solved  by  making  the  substitution  r\  =  cos  0  and  noting  that 
^^2 


sin0 


d^T] 


=  -sin0 


=  -cos0  =-q 


d0  dq  d0 


(A12) 


(A13) 

(A14) 

(A15) 

(A16) 


=  -sin0 


dq 


d^Y 

d0^  d0l^dqd0^ 


d^yrdq^^  ^  d^q  dY 
dq  ^  I  d0  J  d0  ^ 


-T1 


dq 


(A17) 

(A18) 

(A19) 

(A20) 
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Substitution  into  (A12)  gives  Legendre's  equation 

{1-ti")^-2ti^  +  /(/  +  1)L  =  0  (A21) 

dy\  dr[ 

The  only  solutions  which  are  finite  across  the  range  T]  =  [-1,1]  occur  for  /  =  0, 1, 2, .... 
These  are  the  Legendre  polynomials 

Y^P,{t])  =  P,{cosB).  (A22) 

The  left  hand  side  of  equation  (All)  can  be  expressed  as 

2r  —  +  r^^  +  \kl/-l{l+l)]R  =  0  (A23) 

dr  dr  ' 

which  after  substitution  of  ^  =  k^^^r  and  using 


dr  ~  d^  dr  ~  dC 

d^R_  d  fdKdCVC 

“dcUc  dr  ]dr 


-.2  ^ 
■  ^''dC' 


(A24) 

(A25) 

(A26) 

(A27) 


can  be  rewritten  as 

^  +  ~  +  p-'('+l)]|-  =  0  (A28) 

for  which  the  possible  outward  going  wave  solutions  are  j,(0,  or  h,(Q  (Morse 

and  Feshbach,  1953,  p.  1465;  Morse  and  Ingard,  1968,  p.  337).  These  are  the  spherical 
Bessel,  Neumann  and  Hankel  functions  of  order  I,  respectively.  The  spherical  Hankel 
function  is  related  to  the  other  two  functions  via 

h,  =j;+in;.  (A29) 

In  the  internal  fluid  the  range,  and  hence  can  fall  to  zero,  at  which  stage  n;(^)  and 
h,(^)  are  unbounded.  Thus  the  only  valid  solution  is  j,(C)-  On  the  other  hand,  in  the 
external  fluid  where  the  range,  and  hence  can  approach  infinity,  h,(^)  is  a  valid 
solution.  Consequently,  in  the  two  media 
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Appendix  B 


To  solve  for  A  in  equation  (22)  it  is  first  noted  that 

a  =  (o,o,a). 


(Bl) 


Hence, 


VxA= 


1  d 
rsinQdQ 


(a  sine) 


_ 1 

'll 

1 

\--M] 

J  1 

Lrar'  *'J 

(B2) 


and 

VxVxA= 


a 

sine 3(1) 


-('•A.) 


1  a 

rH - — 

r  sine  3(1) 


—  ^ - ^(a.  sin  e) 

rsineae^  *  ' 


dr 


ae 


1  a 
"sine  ae 


(Arsine) 


4 


(B3) 


where  r ,  6  and  $  are  the  unit  vectors  in  the  spherical  coordinate  system.  Substitution 
of  equation  (B3)  into  equation  (22)  gives 


1  a'A^  _  1  3' 
rdr^ 


cj  dt 


("^)  + 


lA 
r'  ae 


1  a 
sine  ae 


(a  sine) 


(B4) 


Solution  of  equation  (B4)  is  again  achieved  through  a  separation  of  variables, 

A^aA^e,^).^?)  (B5) 


which  gives 


1  d^T  1  a' 


c^T  dt^  CL^fdr' 


A.) 


+ 


1  a 
a^r'  ae 


sine  ae 


(a*  sine) 


(B6) 


Both  sides  of  equation  (B6)  are  independent  and  are  set  equal  to  kj.  Thus  the  left 
hand  side  becomes 


d^T  .,2 


dr 


+  kjCjT  =  0 


(B7) 
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the  solution  of  which  is 


The  right  hand  side  of  equation  (B6)  becomes 

1  a  r  2  ^  1  a  r  ]  a  /  .^0,2 

r^arV  Br  r' 30  sm03e'“*"‘^^ 


which  is  solved  by  a  further  separation  of  variables  into 


a,=X{r).Z{Q) 


to  give 


iar2  3X^  22  ^ 

—  —  r—  +kjr^ - — (ZsinG)  . 

Xarl  dr)  Z  aO  sinGaO 


Again  both  sides  of  equation  (Bll)  are  independent  and  are  set  equal  to  /(/  + 1).  From 
the  left  hand  side. 


— fr’A'|+^!,2x  =  /(/+l)X 

dr\  dr ) 


for  which  the  solution  is 


X  =  Gij,{kTr)  +  Hini{kTr) 


since  applies  to  medium  2,  which  has  bounds  not  equal  to  0  or  oo.  The  right  hand 
side  of  equation  (Bll)  gives,  after  some  simplification. 


_  d^Z  „dZ  f,,,  .  1  1„ 

dG'  dG  I  sin'  G 


for  which  the  solution  is 


Z  =  — P,(cos9). 


Hence, 


"  S  ^  P/  (cos G,  )[GJ,  )  +  H,n,{krr,. )] 
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Appendix  C 

The  coefficients  required  to  calculate  C,  and  D^  from  equations  (40)  and  (41)  are 


Yu  =  7^(^Z.,2^r  j/fc,2^)-(^L,2^)'jr(^Z.,2^) 
l-2a 

(Cl) 

Yi2=i  t  (Kab)  ^l[K2b)-(K2b) 
l-2a 

(C2) 

(C3) 

Yi4  =  l{l+l)n,{k-^b)-  l{l+\)kj-bn',{kjb) 

(C4) 

v,5  =%^h;Uu'>) 

2p2 

(C5) 

=  itiKib) 

2p2 

(C6) 

Y2I  =,  ^  {^L,2«)  j/(^L,2«)-(^L,2«)  j;(^i,2a) 
l-2a 

(C7) 

Y22=,  ^  Uz„2«)  n,(/Cz.,2«)-Ut,2«)  n/UL,2«) 
l-2a 

(C8) 

Y23  =/(/+l)j,(v)-/(/+l)/:7-aj;(/:j-fl) 

(C9) 

Y24  =  /(/+l)n,(/:2-a)-/(/+l)/:7-fln^(^2.fl) 

(CIO) 

(M)'  .  ) 

Y26  ”  „  J;l^L,3^/ 

2p2 

(Cll) 

Y31  =  ^L,2^j|(^r.2^) 

(C12) 

Y32  ~ 

(C13) 

Y33  =/(/+l)j,(M) 

(C14) 

Y34  =/(/+l)n,(/:j-/2) 

(CIS) 
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hr  (*,.,(-) 

Pi 

(C16) 

83  =  ' 

Pi 

(C17) 

Y4I 

(C18) 

Y42  = 

(C19) 

Y43  ~ 

(C20) 

Y44  =/(/+l)n;(v) 

(C21) 

Y46  =  ’  j/(^t,3^) 

P3 

(C22) 

(C23) 

Y52  =2^^,2&n;(/:^2^)-2n,(^^2^) 

(C24) 

Y53  =(M)jr(M)  +  [/(/+l)-2]j;(^r^) 

(C25) 

Y54  =  {k-j-b)n"{kj-b)  +  [l{l+l)-2]n,{kjb) 

(C26) 

Y6i  =2/:i,2aj;(^2.,2a)-2j,(^^2«) 

(C27) 

Y62  =  2^2.,2«n;(/:i,2fl)-2n,(/:2.2u) 

(C28) 

Y63  =  (M)jr(M)+W^+i)-2]j,(M) 

(C29) 

Y64  ={*2.a)n"(^2-«)  +  [^(^+l)-2]n;(^2-<^) 

(C30) 

The  spherical  Bessel,  Neumann  and  Legendre  functions  and  their  derivatives  needed 
to  evaluate  the  above  equations  are  listed  in  Appendix  D. 
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Appendix  D 


The  spherical  Bessel,  Neumann  and  Legendre  functions  and  their  derivatives  are 
commonly  computed  using  the  following  recurrence  relations: 


j/(z)  = 


sinz 

z 

sin  z  cos  z 


21-1  . 


for  /  =  0 
for  /  =  1 
for  /  >  1 


cosz  sinz 
I  Z2~ 


lU)  =  Wcosz+  - 


for  /  =  0 

|)sin. 

for  /  =  1 

z 

for  /  >  1 

2 

— yCOSZ 

iriz)= 


6  3  V  6 

- :r  SinZ+ - T 


z^  z^r  ^\z  z^ 


z  z  z  _ 


for  1=0 

for  l=l 

for  />1 


n;(z)  =  i- 


sinz  cosz 
2 


z  z 
2/-1 


n,_,(z)-n,_2(z) 


for  /  =  0 
for  /  =  1 
for  /  >  1 


sinz  ,  cosz 

I  ^ 

z  z 


n'(z)  = 


-7Sinz+  -T —  cosz 
z  \z  zJ 


for  1  =  0 
for  /  =  1 

for  /  >  1 
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- j  cosz— 


- ^  sinz+ 


for  1=0 


for  /=1 


for  />1 


P,U)  = 


^zPm(z)-“P/-2U) 


for  1=0 
for  l=\ 


for  /  >  1 


For  fixed  z  and  varying  I  the  behaviour  of  the  spherical  Bessel  and  Neumann 
functions  is  different  in  the  ranges  0<  I  <  z  and  z<  I  <°°  ■  The  functions  oscillate 
with  slowly  varying  amplitude  in  the  first  range.  In  the  second  range  the  spherical 
Bessel  function  tends  rapidly  to  zero,  whereas  the  spherical  Neumann  function  tends 
rapidly  to  -oo.  In  this  latter  range  the  recursion  relation  of  equation  (D4)  is  stable,  but 
that  of  equation  (Dl)  is  unstable  and  leads  to  increasing  errors  with  increasing  order. 
A  further  difficulty  arises  in  calculating  sin(z)  or  cos(z)  when  z  approaches  a  zero  of 
either  trigonometric  function.  The  rounding  error  associated  in  representing  %  in 
double  precision  is  reflected  in  the  result.  The  error  may  be  amplified  when  z  is  both 
near  a  zero  and  large,  as  an  imprecise  value  for  K  may  be  used  in  the  argument 
reduction  method  needed  for  the  sine  and  cosine  routines.  The  procedures  used  to 
overcome  these  difficulties  are  outlined  below. 

The  spherical  Bessel  function  was  calculated  by  one  of  three  different  methods, 
depending  upon  argument  z  and  order  /.  For  z<10^,  the  ascending  series 
(Abramowitz  and  Stegun  1965,  p.  437) 


j(^)= _ 5: _ _ 

l-3-5...(2/  +  l)  l!(2(  +  3)  2!(2/  +  3)(2/  +  5) 


was  used. 


For  z  >  25  and  /  <  z  + 1 ,  an  asymptotic  expansion  (Abramowitz  and  Stegun  1965,  p. 
437)  was  used  to  obtain  the  first  two  terms.  The  forward  recursion  formula  of  equation 
(Dl)  provided  sufficient  accuracy  for  higher  order  terms; 
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j/(z)  = 


i{'’[/+|.2)sin[c-^]-e[/  +  iz)cos[c-y]}  for/  =  0,l 


2/-1 . 


for  /  >  1 


where 


Z  {-ir[i+-,2m+\pzy 


,  1  ^  (/  +  m)! 

/  +— ,m  = - ; - - 

2  j  m\r{l-k  +  l) 


(DIO) 


(Dll) 


(D12) 


and  ^  is  a  reduced  sinusoidal  argument  related  to  z  by 


r  =2-  — xroundf— 1-0.001935307 179586476925286767  xroundf  — I  (D13) 

^  32  1271  j  \2tzJ 


For  all  other  arguments  and  orders  a  method  primarily  attributed  to  J.C.P.  Miller  by 
Abramowitz  and  Stegun  (1965,  p.  452)  was  used; 


for <1<L 


(D14) 


where  l  is  a  proportionality  factor  found  by  either 


Jl.  (2) 


(D15) 


when  j;  (z)  is  known,  such  as  when  the  recursion  is  continued  down  to  =  0 ,  or  from 


(D16) 


X(2/  +  l)7f(z) 


In  practice  the  summation  in  equation  (D16)  is  continued  up  to  .  The  Ji{z)  terms 
are  found  for  decreasing  I  via  the  recurrence  relation 
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•/,U)=  1 


for/  =  U,+l 
for  I  = 

for/„,  </</-! 


(D17) 


Cax  rieeds  to  be  chosen  with  some  care.  By  trial  and  error  it  was  found  that  provided 
L„  satisfies 


2.1769+0.175271ogz+2,794!xl0-^(log’)-+6.4446xl0'’(logj)^+1.7055xl0'’(logz)‘' 


(D18) 


the  backward  recursion  is  stable  down  to  =  0.  For  larger  values  of  the  recursion 
is  unstable,  was  usually  obtained  from  equation  (D18),  which  was  always  greater 
than  the  maximum  /  needed  to  evaluate  equation  (12).  More  computationally  efficient 
methods  are  outlined  by  Gautschi  (1967)  and  Olver  and  Sookne  (1972). 

Although  equations  (D2)  and  (D3)  for  the  first  and  second  derivatives  of  the  spherical 
Bessel  function  do  not  diverge  as  quickly  as  equation  (Dl),  inaccuracies  in  their 
evaluation  do  nevertheless  occur  when  they  are  used,  especially  as  the  order  increases. 
Better  results  were  achieved  by  using 


2  Z 


(D19) 


jniU)-j/-l(^)  ^  j;-2(^)-2j;U)  +  j;^2  (2) 


(D20) 


The  recursion  formula  of  equation  (D4)  for  the  Neumann  functions  was  used  as  it  was 
stable  for  the  arguments  and  orders  encountered.  To  speed  computations 


2  z 


(D21) 


A  n,  b)  +  ~  ^  n,.3  (z)  -  2n,b)  +  n,.,  (z) 

'  4z'  2z  4 


were  used  in  place  of  equations  (D5)  and  (D6). 
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Appendix  E 


The  integral  over  time  in  equation  (53)  can  be  evaluated  as 

* 


=  +  +24/, -/>]<;<  (El) 


I  Al  [  |’llA((/-/)qT  (  2A(  f(/-/o)  (  ,) 


.  >tA((/.-/ilc[  I  2A(  f) 

+sin  — 7 - X —  S  YZ TT  : - ^  +  Uc~J 

(/.-/o)  n^-/o)l  A' 


.)  Ar  (  |'nAl(/-/i]J  I  2A<  f(/-/o)  (,  A 


.  ItA((/-/ilJ  I  2At  f(/-/o)  /,  ,1 

(/,-/o)  VU-/o)l  A' 


I  A(  [-ixAllX-zilfT  I  2A(  fU-Zo)  /,  -I 


4 


I  At  [-i7tA((/-/f  ][l+i  /Vrt,.  |  2Al  f(/-/o),  ,  ,  , 
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where  C  and  S  are  the  Fresnel  integrals 

and 

and  erf  is  the  error  function 


erf  z  = 


(E5) 


(E6) 


(E7) 


81 


Calculations  of  the  Sound  Scattering  from  Fluid-Filled  Spherical 
Shell  Sonar  Targets 

Mark  Readhead 

(DSTO-RR-0020) 

DISTRIBUTION  LIST 


Director,  AMRL 

Chief,  Maritime  Operations  Division 
B.F.  Wild 

Mark  Readhead  25  copies 

Library,  AMRL  Maribyrnong 
Library,  AMRL  Fishermens  Bend 

Chief  Defence  Scientist  (for  CDS,  FASSP,  ASSCM)  1  copy 

Director,  ESRL 

Head,  Information  Centre,  Defence  Intelligence  Organisation 
OIC  Technical  Reports  Centre,  Defence  Central  Library 

Officer  in  Charge,  Document  Exchange  Centre  8  copies 

Air  Force  Scientific  Adviser,  Russell  Offices 

Scientific  Adviser  -  Policy  and  Command 

Senior  Librarian,  Main  Library  DSTOS 

Senior  Defence  Scientific  Adviser 

Librarian  -  AMRL  Sydney  2  copies 

Librarian,  DSD,  Kingston  ACT 

Librarian,  Australian  Defence  Force  Academy 

Serials  Section  (M  List),  Deakin  University  Library,  Deakin  University,  Geelong  3217 
NAPOC  QWG  Engineer  NBCD  c/-  DENGRS-A,  HQ  Engineer  Centre,  Liverpool 
Military  Area,  NSW  2174 

ABCA,  Russell  Offices,  Canberra  ACT  2600  4  copies 

Head  of  Staff,  British  Defence  Research  and  Supply  Staff  (Australia) 

NASA  Senior  Scientific  Representative  in  Australia 

INSPEC:  Acquisitions  Section  Institution  of  Electrical  Engineers 

Head  Librarian,  Australian  Nuclear  Science  and  Technology  Organisation 

Senior  Librarian,  Hargrave  Library,  Monash  University 

Library  -  Exchange  Desk,  National  Institute  of  Standards  and  Technology,  US 

Acquisition  Unit  (DSC-EO/GO),  British  Library,  Boston  Spa,  Wetherby  Yorkshire  LS23  7BQ,  England 
Library,  Chemical  Abstracts  Reference  Service 
Engineering  Societies  Library,  US 

Documents  Librarian,  The  Center  for  Research  Libraries,  US 
Army  Scientific  Adviser,  Russell  Offices  -  data  sheet  only 
Navy  Scientific  Adviser 

DASD,  APW2-1-OA2,  Anzac  Park  West,  Canberra  ACT 
Director  General  Force  Development  (Land)  -  data  sheet  only 
SO  (Science),  HQ  1  Division,  Milpo,  Enoggera,  Qld  4057  -  data  sheet  only 
Counsellor,  Defence  Science,  Embassy  of  Australia  -  data  sheet  only 
Counsellor,  Defence  Science,  Australian  High  Commission  -  data  sheet  only 
Scientific  Adviser  to  DSTC  Malaysia,  c/-  Defence  Adviser  -  data  sheet  only 
Scientific  Adviser  to  MRDC  Thailand,  c/-  Defence  Attache  -  data  sheet  only 


(DSTO-RR-0020) 
DISTRIBUTION  LIST  (Contd) 


Mr  Jim  Thompson,  MOD  Sydney 
Dr  Stuart  Anstee,  MOD  Sydney 
Dr  David  Blair,  MOD  Sydney 
Dr  Ian  Jones,  MOD  Sydney 
Dr  Ranjit  Thuraisingham,  MOD  Sydney 
Dr  Don  Richardson,  MOD  Melbourne 
Dr  Kurt  Benke,  MOD  Melbourne 
Dr  Roger  Neill,  MOD  Melbourne 

Mr  Paris  Constantine,  Thomson  Sintra  Pacific,  18/112  Talavera  Rd,  North  Ryde  NSW  2113 

Dr  Andrew  Madry,  Thomson  Sintra  Pacific,  18/112  Talavera  Rd,  North  Ryde  NSW  2113 

Mr  Ashley  Deacon,  GEC  Marconi  Systems,  Faraday  Park  Railway  Rd,  Meadowbank  NSW  2114 

DNCSE,  CPl-6-25,  Campbell  Park  Offices,  Canberra  ACT  2600 

MHIPD,  Dept,  of  Defence,  GP2-3-10,  Campbell  Park  Offices,  Canberra  ACT  2600 

MHCPD,  Dept,  of  Defence,  CP2-3-11,  Campbell  Park  Offices,  Canberra  ACT  2600 

MWSCPD,  Dept,  of  Defence,  CP2-2-19,  Campbell  Park  Offices,  Canberra  ACT  2600 

Commanding  Officer,  HMAS  Waterhen,  Halls  Head  Rd,  Waverton  NSW  2060  2  copies 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  UNCLASSIFIED 


REPORT  NO. 

DSTO-RR-0020 

ARNO. 

AR-008-995 

REPORT  SECURITY  CLASSIHCATION 
UNCLASSIFIED 

TITLE 

Calculations  of  the  sound  scattering  from  fluid-filled  spherical  shell  sonar  targets 


AUTHOR(S) 

Mark  Readhead 

CORPORATE  AUTHOR 

DSTO  Aeronautical  and  Maritime  Research  Laboratory 
PO  Box  4331 

Melbourne  Victoria  3001 

REPORT  DATE 

TASK  NO. 

SPONSOR 

January  1995 

M40171 

MHCPD 

FILE  NO. 

REFERENCES 

PAGES 

510/207/0104 

21 

88 

CLASSIFICATION/LIMITATION  REVIEW  DATE 

CLASSIFICATION/ RELEASE  AUTHORITY 

Chief,  Maritime  Operations  Division 

SECONDARY  DISTRIBUTION 


Approved  for  public  release 


ANNOUNCEMENT 


Announcement  of  this  report  is  unlimited 


KEYWORDS 

Sound  Scattering  Binary  Chemical  Agents  Sonar  Targets 

Sonar  Spherical  Shell  Sound  Pressure 

Calibrating  Sonar  Transmission  (Pulsed)  Acoustic  Measurement 


ABSTRACT 


The  theory  of  the  sound  pressure  scattered  from  a  fluid-filled  spherical  shell  immersed  in  a  second 
fluid  is  developed.  An  arbitrary  combination  of  fluids  may  be  used  to  fill  the  shell,  and  both  the  sound 
source  and  the  receiver  may  be  in  the  near  field.  The  theory  is  used  to  calculate  the  pressure  field  in 
both  the  interior  and  exterior  fluids  and  is  extended  to  determine  the  time-dependence  of  the  reflected 
signal  and  the  target  strength.  The  dependence  of  these  quantities  as  a  function  of  frequency  and  pulse 
type;  shell  material,  diameter  and  wall  thickness;  fluid  combination;  temperature;  and  distance  of  the 
receiver  from  the  sphere  is  examined,  and  is  shown  to  be  sensitive  to  most  of  these  parameters.  The 
performance  of  spheres  filled  with  a  Freon-113™  -  ethanol  or  Fluorinert  FC-72™  -  n-hexane  mixture  is 
considered. 
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